Treatment of electron viscosity in quantum conductance 
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PACS numbers: 

In a recent paper Sai et al. [1] identified a correction R dyn 
to the DC conductance of nanoscale junctions arising from 
dynamical exchange-correlation (XC) effects within time- 
dependent density functional theory. This quantity contributes 
to the total resistance through R — R s + R dyn where R s is the 
resistance evaluated in the absence of dynamical XC effects. 
In this Comment we show that the numerical estimation of 
R dyn in example systems of the type they considered should 
be considerably reduced, once a more appropriate form for 
the shear electron viscosity Tj is used. 

Sai et al. 's expression for R dy ", based on electron-liquid the- 
ory [2], is a one-dimensional integral between the two elec- 
trodes 
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where A c is the cross-sectional area, 77 is the shear viscosity 
of the electron liquid, and n is the electron density. The ex- 
ample system we have considered is the metal-vacuum-metal 
(MVM) junction that can be formed by two jellium surfaces 
separated by a distance d [3]. Since this system is translation- 
ally invariant parallel to the surface, the electron density n(z) 
is defined unambiguously. We have chosen rs = 3 for the jel- 
lium electrodes, to allow comparison with the gold-electrode 
systems presented in Ref. 1 . The density is calculated within 
the LDA. 

Sai et al. used a constant shear viscosity r\ c = 
h(k,F /itao) 3 ' 2 1 120 corresponding to formula (4.7) of Ref. 2, 
appropriate only in the high-density weakly inhomogeneous 
limit, and based on the bulk electrode density. Here the quan- 
tities kp = (37r 2 n)'/ 3 and oq are the Fermi wavevector of the 
bulk electrodes and the Bohr radius respectively. We make 
two changes to this. First, we use Formula (4.10) of Ref. 2, 

TJv ~ «/(60r s 3/2 + 80r 5 1 - 40r s 2/3 + 62r s 1/3 ), more appro- 
priate for realistic densities (r s = 0...20 a.u.), [2] in recogni- 
tion of the fact that rs = 3 is not a high density. This alone 
reduces the dynamical resistance by a factor of 5.37 for a den- 
sity corresponding to rs = 3. Second, in Eq. (1) we evaluate 
this viscosity at the local density rather than taking the bulk 
viscosity to apply outside the bulk region. This further reduces 
the dynamical resistance by a ^-dependent factor of 1 . 16-18.4 
(see Table I), particularly for larger d when the dominant con- 
tribution to the integral in Eq. (1) comes from the low-density 
region in the vacuum (see Fig. 1). Thus a more appropriate 



choice of the shear electron viscosity 77 reduces the dynamical 
resistance by a factor between 6 and 98, causing it to become 

d R s R d c yn R? yn 



117 
241 
712 



1.71 
48.5 
530 



0.275 
5.01 
29.4 
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TABLE I: Non-interacting and dynamical resistances per unit area 
(a.u.) as a function of the separation d (a.u.). The interpolated local 
formulation of the viscosity, rj,,, considerably reduces the dynamical 
resistance (R d , yn ) relative to the high-density bulk formulation of Ref. 

1 (R d / n ). 
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FIG. 1: (Color-online) MVM junction separated at d = 3 a.u. The 
electron charge density and the background positive charge density 
are represented by thin discontinuous lines. Important contributions 
to the integral of Eq. (1) (green thick dashed line) arise within the 
vacuum junction, where the electronic viscosity (solid line) is lower, 
tending to reduce the dynamical resistance. The value of viscosity 
used in Ref. 1 is indicated by tj c . 



very small compared with the single-particle resistance in all 
cases studied. 
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